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Abstract 

It has been observed that simultaneous explanation of the solar and at- 
mospheric neutrino deficits and the reported evidence for v„ — ► z/ e oscillation 
from the Los Alamos Liquid Scintillator Detector (LSND) requires at least 
one extra neutrino species in addition to the three known ones. The extra 
neutrino must be sterile with respect to the known weak interactions. We 
present a new mass matrix for these four neutrinos in which the LSND effect 
and the atmospheric neutrino deficit are governed by only one parameter. We 
investigate the phenomenological implications of such a mass matrix ansatz 
and suggest possible ways to understand it in gauge theories. 



I. Introduction 

These are very exciting times in neutrino physics. For the first time in its 
history, there are several hints for a nonvanishing mass for at least two of three 
known neutrinos which look very promising and credible. They come (i) from the 
observation of the solar neutrinos in various experiments and their disagreement 
with the predictions of the standard solar model]!]]: the earlier experiments from 
Homestake, Kamiokande, SAGE and GALLEX0 and the most recent high statistics 
confirmation of these results by the super-Kamiokande experiment [§] and (ii) from 
the observation of the atmospheric neutrinos by several previous experiments^. || 
and the most recent confirmation of the earlier results by the super-Kamiokande 0] 
collaboration. Then there is the result from the Los Alamos Liquid Scintillator Neu- 
trino Detector (LSND) which gives the first laboratory evidence for the oscillation 
of both z/ M — > z/ e || as well as z/ M — > v e typej7j. The existence of nonvanishing neu- 
trino masses can be inferred from these data once one tries to understand them in 
terms of neutrino oscillations which depend on the mass difference square between 
the neutrino species (i.e. 5rnf,-). 

It has been pointed out that since the above processes require three different 
5m 2 , it is impossible to understand them within the conventional three neutrino 
(z/ e , z/ M , v T ) framework; rather, one has to include a sterile neutrino which is ul- 



tralight (denoted by i>s)|| ||. To see this in detail, note that the MSW|T0[ type 
solution to the solar neutrino puzzle including the latest super-Kamiokande results 
yields the following ranges (the so-called small angle solution) for the neutrino mass 
differences and mixing [|TTJ, 



Sm 2 ei ~4x 1(T 6 - 1.2 x l(T 5 eV 2 , sin 2 2fl ei ~3x 1(T 3 - 1.5 x 1(T 2 (1) 

Next, a best fit to the presently available atmospheric neutrino data from Kamiokande, 
1MB and Soudan as well as super-Kamiokande assuming — v T oscillation seems 



to requirepfl at 90% CL: 



4 x 1(T 4 eV 2 < 5m 2 atm < 5 x KT 3 eV 2 , sin 2 26 atm ~ .76 - 1 (2) 



Although the v e — interpretation of this deficit for the older data from Kamiokande 
is ruled out by the recent CHOOZ[|T3|] results, inclusion of the super-Kamiokande 
data still leaves an allowed domain for 5m 2 epu . In this paper, we will assume — v T 
oscillation to be the explanation of the atmospheric neutrino puzzle. Finally, the 
LSND data along with the constraints from the reactor experiment at Bugey|l4] 



and E776 at BNLjlSj suggest that: 



.2 eV 2 < 5m 2 efl < 3 eV 2 , 2 x KT 3 < sin 2 26 e ^ < 4 x KT 2 



(3) 



We thus see that three different 5m 2 values are needed and hence the obvious need 
for an extra sterile neutrino species. 

Once one accepts the idea of four-way neutrino oscillation as a necessity to solve 
the neutrino puzzles, the next question is what is the nature of the mass matrix that 
can accommodate the known data and finally what it implies for physics beyond the 
standard model. It is the goal of this work to discuss the first question followed by 
an attempt to study the resulting theoretical implications. 
II. The four neutrino mass matrix 

A simple mass matrix ansatz that embodies the above mixing pattern was sug- 
gested in[|l^, [17j and analyzed in detail in |L7j . This is given in the basis (u s , u e , v T ) 



by: 



M 



/ /ii /i 3 \ 

fi 3 a 

5m 

\ a m 5 J 



(4) 



where one requires that /13 fi\, 5 a m. A characteristic of this mass matrix 
(to be called henceforth MI) is that the atmospheric neutrino oscillation is governed 



by the parameter 6 whereas the mixing angle in the v e — oscillation is given by 
the parameter a. They are unrelated. 

Our proposal in this paper is to choose instead the following mass matrix: 



M 



(lH H \ 

/i 3 a 

5 m 

\ a m —5 j 



(5) 



This mass matrix is the same as the one in eq. (§) except for the negative sign 
in the v T — v T entry. We will call this MIL We will see that this simple change of 
sign has profound phenomenological consequences. In particular, it implies that the 
same parameter a is responsible for nonvanishing effects both for LSND as well as 
atmospheric neutrinos. 

In order to study the phenomenological implications of this mass matrix, let us 
obtain the eigenvalues and the mixing matrix. We will assume that the above mass 
matrices are defined in a basis in which the charged leptons are mass eigenstates. 

We choose the following hierarchy among the parameters: /i3 <C fi\ <C a <C 
5, m. Note already a difference from the mass texture MI that 5 is much larger. To 
obtain the two largest eigenvalues, we ignore the small parameters ^3 and a. It 
is then clear that the two eigenvalues are equal in magnitude and opposite in sign 



i.e. A 3j4 = ±\/m 2 + 5 2 . This, of course, is not of physical interest since it does not 
lead to atmospheric neutrino oscillation. However, it gives us an idea of how big the 
parameters m and S should be. For instance if z/^ jT are to constitute the hot dark 



matter of the universe, then \/ m 2 + S 2 should be in the eV range. 

To get the mass splitting between the and v T (mass eigenstates denoted by 
z/3 and 2/4 respectively), we clearly have to keep the a and consider the 3x3 matrix. 
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Diagonalizing it, we find that the masses of are now split with 

a 2 5 

m\-m\ = Sm 2 atm ~ ==== = (6) 

V o + m z + a z 

Furthermore, the LSND u e — mixing angle is given (see eq. (|TT|)) by: 



We thus see that in Mil, the single parameter a is responsible for the atmospheric 
mass splitting as well as the LSND v e — mixing angle. This is different from 
what happens in the MI case. We can therefore conclude that if v e — oscillation 
becomes very small and if the 5m 2 inferred from the atmospheric neutrino deficit 
remains where it is, Mil will be ruled out. In fact, an even stronger statement can 
be made i.e. if one takes the 5m 2 values from the LSND experiment consistent with 
E776 and Bugey (see eq. (|)), one has 0.02 < 9 efl < .1. From eqs. (|6|) and (|7|) 
8m 2 atm ~ 9 2 e ^m8( m2 +2 2 )^ . This would at most be 1.3 x 10~ 3 eV 2 , implying that a 
significant part of the allowed Sm 2 tm will be inconsistent with this model. In other 
words, as the atmospheric neutrino data becomes more precise Mil will be testable 
unlike the MI texture where by adjusting the different parameters 5 and a a much 
larger range of values for the above experimental observables can be accommodated. 

Let us now turn to the v e — v s mixing which is responsible for understanding 
the solar neutrino puzzle in MIL For this purpose, we now consider the full 4x4 
mass matrix and find for the u e and v s (mass eigenstates denoted by v-i and v\ 
respectively): 



a 2 5 x 2 



To address the solar neutrino puzzle via the MSW mechanism it is further necessary 

for the sterile neutrino to be heavier than the u e . 

4 



In order to discuss further tests of the model, we write the complete mixing 
matrix among the four neutrinos in the following form: 



U 
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— s 





(9) 



where 



2C ' 

"A*3 



2C ' 



a C-<5 



A* 1 m 2 +<5 2 



a2(5 , C V 2C 



63 - ~ 



a /C + o" 
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c V 2C ' 64 ~ C 2 ' 65 ~ C 2 



(10) 



and C = V m? + 5 2 = ^Sml^. Note that C and c refer to different quantities. 

In terms of these mixings, we can write down the oscillation probabilities for 
various neutrino species P{yi — > Vj) as follows: 



■ v s 

V-r 



~ 1 - 4e 2 sin 2 A s - 4(e 2 + e^) sin 2 A 

~ 4cV sin 2 A 43 

~ 4cse2f3 sin 2 A 43 — 4e 4 (ce2 + ses) sin 2 A 

~ 4e 2 sin 2 A 2 i 



~ 4 

= 



~e 2 e 3 cs sin 2 A 43 + e 5 (-se 2 + ce 3 ) sin 2 A 



with four of the Ay equal, i.e. A 42 = A 



(11) 



41 



where we have defined Ay 

A 32 = A 3 i = A. The remaining two Ay are very different and explicitly included. 

We have used the usual units i.e. 8m 2 is in eV 2 ; L is in kilometers and E is in GeV. 
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The allowed ranges of the parameters where all constraints are satisfied are 
presented in Table 1. 



Parameter 


Allowed Range 


m 


0.39 eV - 1.73 eV 


5 


2.2 x KT 3 eV - 0.85 eV 


a 


2.0 x 10" 2 eV - 0.2 eV 




5.5 x 10" 3 eV - 6.5 x 10" 3 eV 


V>3 


9.0 x 10" 6 eV - 1.3 x 10" 4 eV 



Table 1: The allowed ranges for the parameters of the model consistent with present 
accelerator, reactor, solar and atmospheric neutrino results. 



These ranges must be taken together with the constraints: 

rp' R ()Tn 2 

m = Csm(29 atm ); 5 = Ccos(20 atm ); — = (12) 



uiC 2 — a 2 5 

3.3 x KT 4 < P < 2.1 x 10~ 3 (13) 



The ranges in Table 1 and the above constraints are based on the assumption that 
/ii is positive. If /xi is negative then (a) the constraint fll3|) above does not apply, (b) 
Table 1 refers to the modulus of /ii, and (c) the range for /13 changes to 1.5 x 10~ 4 < 
/i 3 < 7.4 x 10~ 4 . 

To summarize this section, we note that in our mass matrix we use five param- 
eters to accommodate six experimental numbers i.e. three mass difference squares 
and three mixing angles for the solar, atmospheric and LSND experiments subject 
to the constraints from the Bugey and E776 results. (We have checked that the 
results from the CHOOZ and Fermilab E531[18|] experiments are automatically sat- 



isfied when the parameters are in the allowed ranges.) More importantly, because 
of the intimate link between the Q eil and 5m 2 atm) our mass texture Mil is more easily 
testable than the MI texture. 
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III. Other experimental tests 

Turning now to the tests of the model, the first well known point to emphasize 
is that the v e — *> v s oscillation solution to the solar neutrino problem will be tested 
once the SNO measures the neutral current effects of the solar neutrinos. Unlike the 
v e ~^ oscillation scenario, where one would expect &cc/$nc — -4, for the v e — > v s 
case, one should get $cc/$nc — 1- Here $ represents the signal for the appropriate 
process before the incorporation of the cross-sections and detection efficiencies. 

The second well known test of the model is of course in the long baseline 
search for — > v T oscillation in experiments such as MINOS, where the oscillation 
probability is driven by A atm and we expect 



For E — 10 GeV, for the MINOS experiment this would imply a transition proba- 
bility of about 25% for 5m 2 ~5x 10~ 3 eV 2 . 

Let us now turn to v e — > v T oscillation. In this case, there is a difference 
between the MI and the Mil scenarios. In both scenarios there are two mass scales 
that contribute to this probability: the one that governs the atmospheric neutrino 
deficit and the other which is associated with the LSND effect. However, from eqs. 
(H, [11]) it is obvious that in the Mil scenario the contribution of the sin A term is 
considerably larger than in the MI case where 5 <C m. Taking values for 62 — .02, 
£3 ~ —.03, and £5 ~ .03 as dictated by our choice of parameters given earlier, we 
expect that in a long baseline experiment, 



In the MI scenario, the first term is negligible. 

The predicted v e v T oscillations driven by two vastly different mass-scales 

but of comparable strength could be observable in future facilities (presently under 
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(14) 



P(u e -> v T ) ~ 3 x 10- 4 + 10- 3 sin 2 A atm 



(15) 



contemplation) such as the muon collider with a long baseline neutrino experiment [IS 
For the purpose of illustration we consider the detectability of a v T signal through 
the — > fi~ decay where the parent v e beam of E max = 20 GeV originates from 
the decay of /i + s produced for a muon collider set-up at Fermilab. Our results are 
presented in the Figure. We find that a 10 kT detector at Soudan will be able 
to scan, at the single event/year signal level, significant portions of the predicted 
mixing angle - mass- square- difference space of both the mass scales but Gran Sasso 
will be sensitive to a smaller part of only the LSND mass-scale. If a 10 kT detector 
could be set-up at a distance of 10 km from the source then its reach in the mixing 
angle will be much wider and we show in the figure the sensitivity of such a facility 
assuming 100 signal events/year. 

There will also be testable difference between the MI and the Mil scenarios 
in the short baseline experiments. In this case, the A 34 contribution to P[y e — > 
u T ) would be insignifcant. As a result, in the MI scenario, one would expect this 
oscillation to be negligible. On the other hand in the Mil scenario, we get 

P(u e -> u T ) ~ 10" 3 sin 2 A (16) 

This is large enough to be observable with a reasonable r detection efficiency 

In our model, the neutrinoless double beta decay vanishes at the tree level 
since its amplitude is proportional to < m Ue > which is given by 

< m Ue >= EiU^rjim^ = M v&Va (17) 

Here M Vf , Ve is the v e — v e element of the neutrino mass matrix in the weak basis 
and r]i = ±1 is the sign of the i-th eigenvalue. As is clear from eq. (f|), this entry 
vanishes in our case. This also happens in the scenario MI, contrary to the assertion 



in Ref. 17 
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IV. Possible theoretical scenarios 

In this section, we explore possible theoretical scenarios that could lead to the 
mass texture analyzed in this paper. Specifically, we want to reproduce the key 
distinguishing feature of the model which is M V)iV)i = —M VtVt within a reasonable 
scenario that can also roughly account for the rest of the texture. 

Let us consider an extension of the standard model where the fermion sector 
is left untouched but the Higgs sector has several additional scalars and the theory 
is invariant under the extra global symmetry S3 x Z 2 where S3 is the permutation 
symmetry on three objects. The extra Higgs fields included are two triplets (Aj) 
with lepton number L = —2, which acquire a tiny vev due to a Higgs see-saw 
mechanism; four doublets H 12 and with identical quantum numbers to the 
standard model doublet and different quantum numbers under the discrete global 
symmetry group; two iso-singlet charged (77+) fields with L = +2; a doubly charged 
iso-singlet field (k ++ ) with L = — 2 and a neutral L = 2 field a. We will assume 
that < a »< H^ 3A »< A° >~ 1-2 eV. We further assume that < Hf >= 0. 
The S3 x Z 2 transformation property of the Higgs and fermion fields are presented 
in Table 2. 



Fields 


(i) 




('•') 


e c 


H3,4 




(£) 


k ++ 


(:) 


a 


S 3 x Z 2 






















quantum 


(2,+) 


(V) 


(2,+) 


M 


(1, +r) 


(2, +) 


(2,+) 


(1,+) 


(2, -) 


(1,+) 


numbers 























Table 2: The lepton and scalar fields of the model and their x Z 2 quantum 
numbers. There are two one-dimensional representations of S3. Here 1 stands for 
the invariant representation. 
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The Yukawa Lagrangian invariant under the gauge group and the extra global 
symmetries is: 

C = f 1 [2L T L lx A 1 + (L lJl L lx -L T L T )A 2 ] 
+/iiL e # 3 e c + h 2 (L^ c + L t t c )H 3 
+h 3 [(L M r c + L t1 i c )H x + (L^ c - L T r c )H 2 ] 

+f 2 (L T r) 2 + L^)L e + f 3 (r c r c + /i c /i c )^ + / 4 e c e c A;- (18) 

The Higgs Lagrangian will contain two parts: one part consisting only of dimension 
four operators that are invariant under the gauge as well as the global symmetry 
and a part consisting of terms with dimension two or three which will softly break 
the global symmetries. In the equation below, we give only the part of the potential 
relevant to generating the neutrino Majorana masses out of two loop effects. 

£m g9 s = V + fn(r]f + vf)k~~ + A(A 1 // 1 + A 2 H 2 )H 3 a + m 2 A{A 2 

+X'H 4 [(A 1 H 2 + A 2 ifi)77i + (A 1 H 1 - A 2 H 2 )r) 2 ] + h.c. (19) 

Let us now study the tree level mass matrix for the three known neutrinos in 
this model. It is given by 



Mn 



/ \ 

fiv 2 fiVi 
V f\vi -hv 2 j 



(20) 



where Vi =< A° >. Thus this has the desired structure for the scenario Mil in the 
— v T sector. Incidentally, the A term in the Higgs Lagrangian can be used to 
suppress the Aj vev's via a see-saw mechanism. 

The v e — v T entry (i.e. a) in the mass matrix arises via a two loop diagram. 
The magnitude for this term can be estimated to be 

~ (16vr 2 ) 2 M 4 1 ] 
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where M is a typical scalar mass. If we choose /j ~ 4 x 10~ 3 and fi ~ 10M ~ 1 TeV, 
then we get a « .01 eV which is of the right order of magnitude. We have checked 
that the contributions to other entries of the neutrino mass matrix are adequately 
small. It is important to remark that for our mixing considerations to hold, the 
charged leptons must be diagonal. This is guaranteed by < Hi >= 0. We have to 
fine tune hi < H 3 > and h 3 < H 2 > to get the muon and the tau lepton masses 
right. 

So far we have ignored the sterile neutrino. It can be incorporated into the 
model in one of two ways. An elegant way to do this is to double the particle 
spectrum and the gauge forces of the standard model by including a mirror sector 



to the standard model as is done in Ref. [20] and inducing the Majorana masses for 



the neutrinos via non-renormalizable operators. The lightest neutrino of the mirror 
sector will be identified with the sterile neutrino. The non-renormalizable operators 
will in general induce nonzero values in the v e — entry as well as in the v e — v e 
entry but these contributions are in the range of 10~ 6 eV so that they do not affect 
any of our considerations. The mirror v' T masses will arise from triplets whose vevs 
will be larger than the corresponding ones in the visible sector since the see-saw 
mechanism giving their vevs is proportional to doublet vevs, which for cosmological 
reasons must be larger in the mirror sector: < H mirror >~ 30 < H >. 

Another way is to add the sterile neutrino as an extra fermion to the standard 
model and use the radiative mechanism to suppress its mass to the desired level. By 
adding new Higgs fields to the model already described, this can be accomplished. 
We do not elaborate on this here. 

V. Conclusion 

In conclusion, we have presented a new four neutrino mass matrix which fits 
all neutrino observations with five parameters and where the u e — mixing and 
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the Sm^ T are related to each other unlike the simplest four neutrino mass matrix 
discussed in the literature. Thus, in this scenario the atmospheric neutrino deficit 
seen in the present data implies the observable effect in — v e oscillation as found 
at LSND. This model can also be tested in the short as well as long baseline v e — v T 
oscillation searches. We also construct a theoretical scenario which could lead to 
the mass matrix of the type under consideration. 
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FIGURE CAPTION 



The prediction of this model for v e — v T oscillation in an effective two-flavor plot is 
shown as the rectangular regions. The sensitivity of experiments located at Gran 
Sasso (L=9900 km), Soudan (L=732 km) and a new facility (L=fO km) are shown 
when the source is a u e beam obtained from a parent muon beam (see text) at 
Fermilab with E max = 20 GeV. The Gran Sasso and Soudan results correspond 
to a single-event signal/year while for the new facility at 10 km the results shown 
correspond to 100 signal events/year. 
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